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5 ■ Abstract 

T^lj- . In this note I revisit the calculation of partition function of simple one-dimensional systems 

solvable by Bethe Ansatz. Particularly I show that by the precise definition and treatment of the 



QJ ' partition function the nontrivial normalization factor proposed in a recent work to give the correct 

-*— * ' 0(1) corrections to the free energy can be derived in a straightforward manner. 
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I. INTRODUCTION 

Recent developments in both field theory and solid state physics have shown, that in 
certain problems, where surface or impurity effects are important, non-macroscopic contri- 
butions to the thermodynamic quantities like the free energy may play an important role 
[l|, |2|. Among these systems those solvable by Bethe Ansatz (BA) due to the exact treata- 
bility are of special importance. For these one- dimensional (ID) systems the free energy is 
calculated following the method developed by C.N. Yang and C.P. Yang |3|] for the 6 Bose 
gas. The basic idea of this method is that through the density of momenta (rapidities) 
an entropy can be defined and a free energy as a functional of the rapidity density can be 
constructed. The minimization of this functional yields both the equilibrium state and the 
macroscopic part of the free energy of the finite temperature system. Based on this idea 
the present author developed a method to calculate the 0(1) free energy contributions of 
the states near to the equilibrium ^y], however the contributions found have not met the 
expectations J5[: for periodic boundary conditions (PBC) no 0(1) corrections have been 
expected, but the calculation gave some, and also for the case of open ends with integrable 
boundaries (IB) some of the obtained terms were not of the expected structure. Recently, 
based on intuitive arguments concerning the density of states in the configuration space 
a nontrivial normalization of the partition function has been proposed [6|, by which these 
differences can be dissolved: in the case of PBC the corrections are canceled, while for the 
IB case they are corrected. In the present note we derive this nontrivial factor directly by 
the careful definition and treatment of the partition function. This way this work confirms 
the proposal of |6| and completes [^ . 

As a starting point we briefly review the ingredients of the calculations. This will 
serve also to make clear our notations and expose the problem in a more tractable form. It 
has been tempting to formulate our treatment in a general form, as however the derivation 
of the entropy term for the free energy functional is slightly different for closed and open 
boundaries we treat the two cases separately: first we deal with the case of periodic boundary 
condition in more details, and in case of integrable boundaries we point out the differences 
only. 



II. BA AND BA THERMODYNAMICS 



We consider a system with BA equations 

N 

Lp{9,)+J2m-0j) = 2nh. (1) 

Here the 6 are the rapidities of the particles, p{6) is the momentum of a particle with 
rapidity 9, the 0(6* — 9') is a phase closely related (up to a constant equal) to the phase 



shift arising when a particle with rapidity 9 is scattered on an other with 9' , and the / 
quantum numbers are either integers or halfs of odd-integers depending on the number of 
particles A^. (For later purposes we chose the Riemann sheets to have continuous at zero 
argument.) To each set of quantum numbers {/,-, Ij ^ //, if j 7^ O defines a set of real 
rapidities {9j,9j 7^ 9i, if j 7^ /}. The wave functions belonging to the different solutions of 
([T]) are orthogonal and form a complete set. The energy of such a state (modified due to the 
chemical potential /i if necessary) is the sum of the contributions of the individual particles 

N 

E~fiN = E{9,, 92,...,9n) = Y1 <^i)- (2) 

«=i 

The finite temperature description ^] is based on the idea, that for a macroscopic 
system the roots of ([T]) can be described by densities, and the thermodynamic quantities 
can be given by these. In particular the p{9) density of particles is defined so, that the 
# of ^j G {9,9 + /S.9) is Lp{9)A9 and the ph{9) density of holes is given by the equation 
obtained from ([1]) 

p{9)+p,{9) = ^^ + lK{9-9')p{9')d9' with K{9) = ^^. (3) 

All p{9) functions are physical, for which ([3]) yields a non-negative Ph{9), as for these densities 
Jj quantum number sets can be constructed which define 9i roots distributed according to 
p{9) (with a certain accuracy). The number how many ways this can be done, i.e. the number 
of states represented by one single rapidity density is estimated by the combinatorial factor 

Through this the entropy is 

S[pi9)] = \nn[pi9)], (5) 

and a free energy functional can be defined: 

F[p{9)]=E[p{9)]~TS[p{9)] with E[p{9)] =J2e{9)p{9)A9, (6) 

and T being the temperature. The macroscopic part of this is 

FL[pm=Lje{9)p{9)d9- 

L fT[ip{9)+pH{9))\n{p{9)+ph{9))-p{9)\np{9)-pt,{9)\nphmd9. (7) 



The equilibrium density po{9) is given by the minimization of (ITj) under the constraint of 
([3]). This gives that in addition to ([3]) the densities satisfy 



too, where (3 = l/T, and 

e{e) = e{e) -T j K{e- e') in (l + e-^^^^')) dO'. (9) 

Finally the macroscopic free energy is the minimal value of Fl. 

In j^ we calculated corrections to this macroscopic free energy by calculating the sum 

E e-/3^[P(^)] (11) 

all pie) 

with an accuracy enough to have the next to leading order corrections too. This led to 

^ e-/3F[p(^)] ^ g-/3F^,„ + A5^ (^2) 

all p{B) 



but according to Pozsgay's proposal [6[ the correct partition function is of the form 

Z=U Y. e-/^^[/^(^)] leading to z = U e' 1^ ^-^^^ + ^^ (13) 

all p(0) 



with M being a well defined 0(1) factor to be discussed later. The important claim in ( IT3|l 
is, that f lTT]) itself reproduces the partition function up to a normalization factor only. To 
clarify this we recalculate the partition function starting from the very beginning, as we 
think, that in a correct treatment all normalizations are defined uniquely in a "natural" 

way. 



III. PARTITION FUNCTION AT PERIODIC BOUNDARY CONDITIONS 



The grand canonical partition function is 



oo 



Z=Y.Zn (14) 

N=l 



Zn= E exp{-/3E(ei,02,...,M}, (15) 

/l</2<...<-fjV 



with 



and our aim is to transform fll4tll5p into a functional of the rapidity densities used in the 
thermodynamical description. The steps of this are the following. 




1. First we write the sum in (fT5|) in a form, which can be well approximated by an 
integral. The problem in this is that the requirement /j < Ij is essentially different 
for the cases with the //s being discrete or continuous variables (see Appendix Rl). 
thus one has to represent the S/i</2<...</iv i^ ^ form, which does not contain such 
restrictions. This is possible if, as in our case, the summand is symmetric in the 
variables. The idea is, that from an unrestricted sum we subtract the contributions of 
the nonphysical (Jj = Ij type) configurations. The technical details of the procedure 
are given in Appendix El here we give the result only: 

Y: fih,l2,..-lN)=T.<P) E f^''Kh,l2,...In). (16) 

/l</2<..-</iV P h,...,In 

Here the summation on the right-hand side goes over the partitions P of the number 
A^, and in f^^' groups of the Ij parameters corresponding to the elements in P are 
equal: for example if a partition is P : {pi + p2 + . . . + Pn}, then 



(17) 



In the partition P : {1 + 1 + .. . + 1 = A^xl} the summation is taken without restriction 
over the J's, and the coefficient is a(-P) = ]^. For the partition p:{2 + l + ... + l = 
2 + (A^ — 2) X 1} in the summation the first two variables of / are kept equal, and the 
coefficient is a{P) = ]^(2)- ^ systematic way to construct the a{P) coefficients and 
a proposal for their form is given in Appendix [A] 

The summand / depend on the Ij quantum numbers through the 6i variables deter- 
mined by ([T]). For this we have to define the (nonphysical) solutions of ([T]) for I sets 
containing equal quantum numbers too. In order to do this we consider (2.1) as an 
analytic relation between the elements of the sets {Ij} and {Oi}, and for equal //s 
we shall take the solutions for ^'s obtained through the limiting process in which the 
corresponding J's tend to the required value. In practice this means, that equal J's 
will define equal 6''s, and each 6j will be summed according to its multiplicity, i.e. for 
a partition P : {pi + p2 + . . . + p„} 

n 

Lp{ei) + Y. <Pi(^^ - Qi)Pi = 2vrJ.. (18) 

'J 

2. Next we replace the summations by integrals: 

m:---T.^jj---jf{dh. (19) 

-fl h In « = 1 

For this from now on we have to consider the Jj quantum numbers in ([1]) and ( !T8|) as 
continuous variables. 

As in our derivation this integral representation of the sums over the quantum numbers 
has a central role, we have to note on the error introduced this way. It seems plausible, 



that this step is accurate enough, but to give a quantitative estimate on the error needs 
different considerations for the different models. In Appendix [E] we analyze the case 
of the S Bose gas with PBC in more details. 

3. The integration variables are changed from the J's to the ^'s 

n<//.^//.../detg}n.<,, (20) 

The Jacobi determinant is obtained from f llSp : 

det jl^l = det{Lg{e,)S,, - K{e, - e,)p,} (21) 

where 

t j = l 

(with K(9) given in (j3])). Factorizing (12T|) as in [6i] one has 

detl^§-^=(uLg{9,)]D(^\e„...e^) (23) 

with 

It is important to note, that both g{Oj) and det {6ij—K{9i — 9j)pj/LQ{9i)} are contin- 
uous in the 6''s. (Although the dimensions of the matrices {Sij — K{6i — 9j)pj/Lg{9i)} 
are different for the cases 9k = 9i and 9k — t- 9i, the corresponding determinants are 
equal.) Due to this, if we denote by D{9i,92, ■ ■ ■ ,9iy) the determinant belonging to 
the case of all Ij (i.e. all 9j) different, then D^^\9i, . . . 9n) can be obtained by taking 
groups of the variables equal like in ( TT71) . This way we have 

E f{Iu...lN) = 

h<l2<-<lN 

E a{P)j. . .//(^H^i, . . . 9^,)D^^\9,, . . . 9n,)\^ LQ{9,)d9A , (25) 

where the a{P) coefficients are the same as in flT6|) . 

4. We divide the 9 axis into A9a intervals. If this intervals are small enough (as re- 
quired according to j^), all /(6'i, . . . , 9^), -0(6*1, • • • , ^n) and (){9) can be given accurate 



enough by giving the Na numbers of ^i's falhng into the interval A6a {Na = # of ^j G 
A6a)- For this we may write 

where 9oc is a mean value in A9a, and 

f{9,, ...,9n)-^ /({iVj) = exp {-(3E{{Na})} , (27) 

E{{N^}) = J2<^c.)Na, (28) 

a 

D(gi, . ■ ■ , g^) ^ Di{N^}) = det 1^., - ^^^2g(ff ^" }• (2^^ 

We note here, that for A'^^ each 9i in A9a have to be taken with the multiplicity pi, 
thus J2a ^a = N. We note also, that as all variables within a single interval A9o, 
are taken equal to the same mean value 9^, distinguishing f{{Na}) {D{{Na})) with 
respect to P has no meaning any more. 

5. The 9 integrals are completed under the restriction, that the number of ^'s in A9a are 
Na- Approximating q{9j) in \Tlj^iLg{9j)d9j) by the same constant Q{9a) H 9j e A9a, 
and following the arguments of Appendix O leading to (1C6P (taking also into account, 
that now due to (l25l) each d9i belonging to the integral over A^^ is multiplied by 
Lq{9q)) the result of the integrals over the intervals A9a is 

thus 

Zn= E exp {-m{N.])] D{{N^}) J] M^^^^A , (31) 

where Y,\y^ n^=n\ means summation over all possible A^^ sets consisting of non- 
negative elements satisfying Y^ N^ = N. Finally 

Z = ...E... exp{-/3E({Arj)}D({Arj)nf^'^^;^^M. (32) 

We emphasize, that the factor ( l30l) appearing in the partition function is only formally 
a combinatorial factor, as it is rather the result of some 9 integrals. (For this for 
example LQ{9a)A9a need not to bee an integer.) 

6. At thispoint we modify the notations to make them coherent with those of the earlier 

works [41]: we define the densities according to 

N^ = Lp{9o,)A9^, (33) 



and introduce the density of holes according to 

Q{eo) = p{e^) + Ph{e^). (34) 

With this we arrive at the (partly) familiar expressions 

p{e) + PH{e) = ^^ + EKi0- 0M0a)Ae^, (35) 

ZTT OU Q 

E{{N^}) =LJ2 e{e^)p{e^)M^ = E[p{e)l (36) 

a 



and 



n[p{e)] = \{u{p{e^),pn{e^)). (39) 



Defining the free energy functional in the usual way (see ([6])) we arrive at 

Z= Y. e-(^^^P^^)^D[p{d)] (40) 

with Z){p(6Iq)} nieaning summation over all the p{9a) sets obtained from the {No} sets 
through ( l33l) . 

Actually this (jlOj) is the expression we wanted to derive. It differs from fITT]) (what 
is practically the starting point of |^) in the factor of D[p{9)]^ meaning that contrary to 
the naive application of Yang's ideas |3[ ( ITTj) is able to give the partition function up to a 
normalization factor only. 



From this point the calculation goes as in [4|. We replace the sum over {piOa)} by 
integrals: 

E ^ f...fl[iLA9^M0a)), (41) 

{pidc)} •^ -^ " 

thus we have 

Z= f...fe'^Pip(^)]D[p{e)]l[{LAe^dp{e^)). (42) 

Next in the free energy functional F the entropy term lnr2[p(^)] is expressed by Stirling's 
formula containing also the terms next to leading order. The macroscopic part of F[p{6)] is 
expanded up to second order in r{9) = p{6) —po{9) around the equilibrium density Pq{9), and 
the Gaussian integral obtained this way is evaluated. In this procedure the sub-macroscopic 

8 



terms of the entropy (that in fact regularize the functional integral) and the determinant 
D[p{9y\ can be taken at Po{0). The latter taken out of the integral yields 

Z = D[pom J... J e-/5^[/^(^)] n (LAe^dpie^)) (43) 



The evaluation of the integral in fH^ is given in details in ^ and is recited also in |6[ , thus 
we do not repeat it here, just cite the result: 

Z = ATe"^^™" + ^^. (44) 

Here F^in is the macroscopic free energy given by flUJ]) . AS* is the correction due to the 
contributions of the states near to the equilibrium (saddle point fluctuations) 

and 

J\f = D[pom (46) 

is actually the same as the normalization factor proposed in |6|: due to ([8]) 

^[,„„„.,et{.„,-«^). (47) 

As a consequence, although the origins of A/" and e^*^ are completely different, for the PBC 
they cancel each other indeed. 



IV. INTEGRABLE BOUNDARIES 



In case of integrable boundaries the system is described by the BA equations 

N 

2Lpi9.^ + MO.) + Vde,) + J2 (<^(^^ - Oj) + HOi + Oj)) - 0(2^.) = 27r/,, (48) 

where the v'o/l(^) are phase shifts due to the elastic scatterings on the ends at resp. L, and 
the li quantum numbers are always positive integers. Now only those solutions are physical, 
in which all (real) 6''s have different modulus, and none of them is zero. Also in this case 
the energy is given by (|2]), but now 

oo 

Z=J2Zr, with 2^= J2 exp 1-/3^(^1, ^2,..., Ml- (49) 

iV=l l</i</2<---</jv 

In calculating the partition function we follow the same program as for the PB case, 
but some steps of the calculation have to be modified, as this case differs from that of 



the PBC in two points, both affecting the 0(1) corrections. One of these is the different 
structure of the BAE, the other is that the summation is restricted to the positive integers 
only. 

Due to the modified structure of the BAE a different determinant will appear in the 
partition function. When allowing for the quantum numbers being equal and treating them 
as continuous variables the equation connecting the Ij and 6i variables are 

n 

2Lp{e,) + ip^{e,) + ipLm + ^ (0(^i - e,) + m + e^))p^ - 0(2^0 = 27r/,. (50) 

i 
This gives a Jacobi determinant 

det jl^l = <^ei{Lm)5^, - K-{e,,e,)p,} (51) 

1 " 

idp{e) 1 dM0) , 1 9^0) 2 



with 



where 



and 

K^{9,9') = K{9-9')±K{9 + 9'). (54) 

Factorizing it in analogy with (123|) leads to 

np 
D(^)(^l,^2,...^np)n^^(%)' (55) 

i=i 
with 

D^^\9„ ...9^) = det l^,, - ^^^1^} • (56) 

This appears in the final expression for Z in the form 

with 

p(e) + p^(^) = a{9) + 5: K+(^, 9o)p{9^)M^. (58) 

The restriction < Ji < /2 . . . affects the combinatorial factor entering into the free 
energy functional, and the careful inspection of it is needed to have the "regular" 0(1) terms 
correctly. Due to this restriction the proper representation of the sums by integrals reads 
(see Appendix [D|) : 

E E--- E ^JJ...Jf[{i-smdi,, (59) 

/l=l /2 = 1 In = l *~^ 

10 



with S{I) being the Dirac 5-function, and it is understood, that its integral on the hmit of 
the integration interval is 1/2. Changing the integration variables from the /'s to the 6''s, 
and taking into account that 6{I) = 6{9)/{Lg{9)) we get 



JJ...Jf[ii-5mdi, 



*-i 



np 



//•••/ n (Lgidj) - m)) d9, D^^\e,, 02,... 9np). (60) 

\J=i / 

Completing the integrals by the A^^ intervals we end up at the combinatorial factor 

mo)] - ("<^<*»' "/;'<*:'>r" - '"] n - (P(o.), pm) (61) 

V Lp{eo)Aeo J ^^Q 

with u{p,ph) given by (138|) and the a = index referring to the A^ interval starting at the 
origin. With this the free energy reads as 

F[p{e)] = E[p{e)]-T\nn[p{0)]. m 

As a result the partition function is of the form 

Z= Y. e-^^[/'(^)]5[p(^)], (63) 

that is in complete analogy with f l4Up and can be treated the same way leading to 



Z = ATe"^^™™ + ^^. (64) 

Here Fmm is the bulk free energy modified by the "regular" 0(1) corrections due to the 
reflections on the boundaries and the exclusion of the zero rapidity: 

Fmin = Fmin + AF + 0o + 0L (65) 

with 

AF = - In ("l + e~^^(0) )+2T f K{2e) In (l + e~^^(^) ) dO, (66) 

o 

and 

^0/. = -^/%^lnfl + e-W«)),9; (67) 



27r J 89 



the AS* correction due to the saddle point fluctuations is given by 



-1 



11 



and (as ([8]) with the same e holds for po(^) and Ph,o{0) too) the factor M is 

Ar = B[p„(9)l=det{«<„-:^-i^^l|^j. (69) 

In this case the normahzation factor Af and the contribution of the saddle point fluc- 
tuations e^*^ do not cancel each other, but after a straightforward manipulation lead to the 
expected boundary independent correction 



A^ + InAf = 

^n/ •••yn7-^k(^i + ^2)^(^2 -^3)... i^(^„-ei) (70) 



" ^ — oo — oo 



=1 Vl + e^ 



indeed. 



V. CONCLUDING REMARKS 

In the present work we revisited the calculation of the partition function for certain 
BA solvable models at periodic boundary conditions and also at integrable reflective bound- 
aries. Our motivation for this was, that the earlier works starting from the free energy 
functional given in terms of rapidity densities, and based on the calculation of the saddle 
point fluctuations [4] turned out to miss important 0(1) corrections to the free energies of 
these systems [6|. To derive the missing pieces we started by the exact expression of the 
partition function given as a sum over the possible quantum numbers entering into the BA 
equations (( ITSll and (H9l)). In several steps (representing these sums by integrals over the 
variables obtained by considering the quantum numbers as continuous variables, changing 
the integration variables to the rapidities, and completing this integrals in a special way) we 
transformed this expression into a one given as a sum over the possible rapidity densities. 
This way we have shown, that the correct expression for the partition function is the one 
given in terms of the rapidity densities in the usual way, but completed by a density depen- 
dent normalization factor arising due to the Jacobi determinant connected to the change 
of integration variables (( H0|) . (l63|) ). This factor (being 0(1) in itself) can be taken at the 
equilibrium density to yield the nontrivial normalization factor proposed intuitively in |6| 

im,m)- 

An interesting feature of the calculation is that the combinatorial factors giving the 
entropy part of the free energy functional are obtained as the results of some integrals over 
the rapidities themselves. As these integrals are slightly different for the PBC and the 
IB case, we have presented the calculation for both cases. The main point, however, the 
appearance of the nontrivial normalization has the same origin, namely the change of the 

12 



variables from the quantum numbers to the rapidities, and inspecting the shght differences 
in the integrals are needed not to miss the "normal" 0(1) corrections in the IB case. 



Since the fundamental work of C.N. Yang and CP. Yang [3| most of the thermodynamic 
descriptions of Bethe Ansatz systems have been formulated in terms of the rapidity densities. 
A different approach of the problem was presented by Kato and Wadati |7|, who calculated 
the iV-particle cluster integrals based on the partition function (I14til5p . and have found 
complete agreement with the results of thermodynamic Bethe Ansatz. It seems, however, 
that such an equivalence exists on the macroscopic level only, and the thermodynamic Bethe 
Ansatz formulated in the usual way is not accurate enough to give the sub macroscopic 
contributions correctly. To be definite, the combinatorial factor (jlj) 

is able to give the number of states represented by a rapidity density p{6) in leading order 
only, and the expression giving also the next to leading order correctly is of the form 



Here the appearance of D[p{6)] results in the normalization factor proposed in [6|. 
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Appendix A 

Our aim is to find an accurate enough integral representation of a sum of the type 
S/i</2<...</iv /(-^i' -^2, • • • -^Af) with symmetric summand /({/j}) and J/s being integers or 
half-odd-integers. The replacement 

N 

J2 fih, /2, . . . In) =^ / /(/i, /2, . . . In) U dh (Al) 

h<l2<...<lN Jll<h<-<lN i=l 

is obviously too rough for our purposes: if we suppose, that the Ij are the half-odd-integers 
between zero and some integer M , and we take f{{Ii}) = 1, then the value of the left-hand 
side is (jy), while the right-hand side (with limits zero and M) gives M^ /N\. The difference 
is caused by the improper treatment of the exclusions Ji ^ I2, I2 ^ I3, ■ ■ ■ In-i 7^ In- Now 

13 



we present a systematic way to solve this problem. As the treatment of an /j < Ij restriction 
in an integral gives different result as in a sum, first we represent the left-hand side of (lAip 
in a form, which does not contain such restrictions, and in a second step we analyze the 
properties of the integral obtainable this way. 

We start with the sum 

So = ^EE---E/ai,/2,.../^) (A2) 

Ii h In 

in which the contribution of the configurations with all Ij being different is the same as 
required, but it contains also the contribution of the terms with some of the /j's equal. For 
this we subtract a sum in which the the configurations with different /'s do not contribute, 
but it eliminates at least the unwanted /,; = Ij contributions: 



^^ ■\'^ / h h In-1 '^~^ ' 



(A3) 



In Sq + Si the configurations with exactly two J's equal do not contribute, but it still counts 
the contributions of the configurations with more than two /'s equal. Next we eliminate from 
So + Si the contributions of the configurations of the types Ii = Ij, Ik = h and Ii = Ij = Ik 
by adding: 



+ 2^ (^) E E ■ ■ • E /(ii:ir[i^^2, . . . /iv-2). 

'V / /i h In--?. 



(A4) 



In Sq + Si + S2 none of the configurations, which are described by one or two equality do 
not contribute. The configurations described by three equalities are of the type i): Ii = Ij, 
Ik = h and Im = In, h): Ii = Ij = h and /; = /,„, or iii): k = Ij = I,^ = I^. These can be 
eliminated by adding 

• • • E /(-^i' -^1' -^2, h, h, h, h, ■ ■ ■ Ins) 

Ins 
-2^(^)(^2~^)EEE---E/(/i>A'V2'^2,/3,.../^-3) 



61^(4) EE- ■■ E fihJuhJul2,...lN-s) 
^^■\^/ h h In-s " 



(A5) 



Now the systematics is clear: in each step we eliminate the contributions with the minimum 
number of equalities. For example in 

n 

E S^ (A6) 



i=0 
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the contribution of the configurations with all /'s different is as required, but it contains 
unwanted contributions from configurations in which groups of J's containing all together 
more than n + 1 elements are equal. The E„_|_i is chosen to eliminate the contributions of 
those configurations, which can be characterized by n + 1 equality signs " =" . The procedure 
is in each step combinatorially well defined, and it leads to the form 

Y: f{h,h....lN)=Y.<P) E /(^H/l,/2,.../np). (A7) 

/l</2<-</iV P h,l2,-,Inp 

Here the summation on the right-hand side (r.h.s.) goes over the partitions P of the number 
A^, the a{P) coefficients are constructed by the above procedure, and in the functions f^^^ 
groups of the /j parameters corresponding to P are equal: for example for a partition 

P ■ {P1+P2 + ■■■+Pn} 

N 



f^'Kh, h,... In) = f{Iu_^,l2^_^, . . . ,In^_^). (A8) 

Pi P2 Pn 

We note here, that the flA7p representation of the sums over the quantum numbers 
is closely related to the one introduced in [7], but it is formulated in a less abstract way. 
Comparing the two yields an expression for the a{P) coefficients: 

aiP) = ^FiP)CiP) , (A9) 

where 

n 

FiP) = l[i-^Y'"''\p^-^y-, (AlO) 

as given in J7|, and C{P) is the combinatorial factor giving the number of ways the sets of 
Pi, P2, ■ ■ -Pn elements can be chosen out of the N ones irrespective of their order. If in the 
partition P : {pi + p2 + . . . + Pn} the element pj is present Uj times, i.e. J2j ^jPj = ^ and 
J2j ^j = n, then 

c(P) = mn^, (All) 

and 

a{P) = n(-l)'^^^"^"'^T77^ • (A12) 



J 



Ujlpj 



In our work we do not use the explicit form of a{P), only the fact, that such a representation 
exists, and some properties of this representation discussed below are exploited. 
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Appendix B 

An important property of the (1A7P representation can be obtained by applying a 
special function for /({/}): 

fIL \ f 1, if /<M, 

f{h,...,IN)=[l[HI^))9iIu...JN) with h{I) = { (Bl) 

\i=i J [ X, if / > M, 

where x is an auxiliary variable and M is a suitable number between the lower and upper 
limits of summations, not equal to any of the possible /'s. Substituting this both sides of 
f lA7|) become polynomials of the variable x: 

1=0 \Ii<l2<---<Ii<M M<Ii+x<h+2<-<lN 

Y^a{p) I Y. +^"^E I ••• f E +^""-E I g^'Kh^h.-.-in,). (B2) 

P \h<M M<hJ \l„p<M M<I„pJ 

Equating the coefficients of the same x powers we obtain relations of the structure 
Y^ Y g{h,l2,...lN)= E a{P)c{P,,P2)x 

i'i<...</jVj<M M</jVj+i<...</jv P=Pi®P2 

E E ^(^^®^^)(/i,/2,... /.,,+„,,). (B3) 

Ii,...,I„p_^ <M M</„p +i,...,7„p _|.„p 

Here Ni is an integer less than A^, Pi and P2 are partitions of the numbers A^i and N2 = 
N — Ni, respectively; P = Pi (B P2 is the partition of A^ emerging as a composition of 
Pi and P2, and J2p=Pi®P2 means summation over all of these partitions; c(Pi , P2) is a 
combinatorial factor giving the number of ways P can be split up into parts Pi and P2; 
finally g^^'^'^^^^Ii, I2, ■ ■ ■ Inp +np ) is the same as g^ ^ just the variables are permuted to have 
the ones corresponding to Pi (those which should be less than M) appearing at the first np^ 
positions. This shows, that if on the l.h.s. of (1A7I) a restriction like in ( JB3I) is imposed, only 
those terms of the r.h.s. contribute which are compatible with this restriction in the above 
sense. This is, however, true on the opposite way around too: if on the r.h.s. of ( ]A7p only 
those terms are taken into account in which the variables Ij counted with the multiplicities 
Pj can be grouped into sets of elements A'^i resp. N2 = N — Ni, and the summation over the 
variables belonging to different sets is carried out over nonoverlapping regions like in the 
r.h.s. of (IB 3 1) , then the result will correspond to a partial sum like that on the l.h.s. of (1B3I1 . 



f lA7|) can be applied to the Ii < . . . < I^i < M and M < Ini+i < ■ ■ ■ < In sums on 
the l.h.s. of (1B3I) . This leads to the relation 

a(Pi)a(P2) = a(Pi © P2)c(Pi, P2), (B4) 
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that can be checked directly for the coefficients given exphcitly in Appendix \^ and it holds 
also for the general formula (IA12p 



Appendix C 



Now we replace in flA7l) the summations by integrals 



np 



/(^)(ji,j2,...,/.,)n^^- (ci) 

il<i2<...<i]V t- i=l 

From the fact, that for certain cases the above replacement is exact, some important prop- 
erties of the integral of the right-hand side can be derived. 

Suppose, that the Ij numbers are half odd-integers between zero and some integer M, 
and /(/i, ^2, • • • , ^n) = 1- The substitution (ICip is exact, if the integral limits are and M. 
From this we have 

As, however, the structure of the polynomial on the l.h.s. is independent of the integration 
limits, 

for any y — x = L. 



In a similar way, we may require that in the l.h.s. of (ICip the first A^i of the /j's should 
be less than the integer Mi, and the rest of them fall between Mi and M (like in Appendix 
iBl) . In this case only those partitions on the r.h.s. contribute, which can be split up into 
partitions of A^i and N2 = N — Ni (i.e. P = Pi® P2 with Pi and P2 being partitions of A^i 
resp. A^2)- The integration limits should be taken as and Mi for the variables belonging 
to Pi, and they should be Mi and M for the others. The result of the integration is a 
polynomial of Mi and M2 = M — Mi which now is 

e;)0 

It follows from this, that if we take the integral such a way, that A^i variable should fall into 
an interval of length Li and the rest into an other nonoverlapping (or even disjoint) one 
with length L2, then the result is 

(30 

This property can be generalized to any number of nonoverlapping intervals: if the integra- 
tion is carried out under the restriction, that A^^ of the variables (counted with the proper 
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multiplicity) fall into the interval Lq,, a = 1, 2, . . ., then the result is 



nfe)- <«) 



Appendix D 

In case of the integrable boundaries we have to calculate sums over the li quantum 
numbers taking integer values with restrictions {0 < /i,/i < li+i}- The restriction of the 
type li < Jj+i can be dissolved as in the general case leading to 

Y: /(/i,/2,.../;v) = E«(^)E E---E f^''\h,h,...K), (Di) 

0</l</2<...</iV P 0</l 0</2 0<I„p 

with P being the partitions of A^. If M is an integer, 

E/(.../...)^ / fi...l...)dl. (D2) 

^=1 1/2 

This formula, however, is not convenient for us, as changing the integration variables to 
the rapidities 6, the limit of integration replacing the 1/2 will be a function of the other 
rapidities. For this we use an other formula 

j:f{...I...)^ f{...l...)dl--f{...l = 0...), (D3) 

7=1 i ^ 

which, by taking into account, that the integral of the Dirac (5-function on the end of the 
integration interval is 1/2, reads 

M ^+V2 

j:f{...I...)^ ff{...l...){l-6il))dl. (D4) 

1=1 i 

This way we have 

M 
E /(/l,/2,.../7v) = 

0</l</2<---</jv 

M+1/2 M+1/2 
E«(^) J ■■■ J f^''\h,h,...In,)Ui^-^i^^))dI^■ (D5) 

For a constant /(.../...) the integrals on the r.h.s. can be completed, and due to l\C2\} they 
yield 

E «(^) / • • • / n ^1 - '^(^^)) ^^^ = f ^ 'n^'^] (^6) 

p ^ i i=i v ^^ / 

for any upper limit L. 
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Appendix E 

The quantitative estimation of the error introduced by the replacement (fT9|) 

h h In *=1 

may need different considerations for different models. In this Appendix we deal with the 
6 Bose gas with PBC, for which the rapidity variables are the wavenumbers (momenta) of 
the particles {9i = ki) determined by the equations (see f lTSj) ) 

Lki + Y,2 tan-i (tLZh\ p^ = 277/,, (E2) 

the energy is quadratic, i.e.: 

e(A;,) = k^-fx (E3) 

(with /i being the chemical potential), thus the f^ function is 

/^(/i, l2...In)= exp I -(5 J2 P3 [k] -li)\. (E4) 

We use the Poisson summation formula (see also 

E 9{I)= E J9{I)e'^'''dI. (E5) 

I=integer J=integer 

According to this, the (fT9l) (or ( lEll) ) type representation of the infinite sum is equivalent to 
taking on the right-hand side the J = term only, and the error introduced can be estimated 
by the neglected J 7^ contributions. For each / we apply flESP (for the sake of simplicity 
we suppose, that the I/s are integers), and we calculate a general term 

„„„|n njn 

... exp 1-/3 J2p, [k] - ^) + 27^1 E hJ, n dh- (E6) 

■'■''' y i=i j=i J i=i 

This, after changing the integration variables from the /'s to the fc's, due to the tan~^ form 
of the phases reads 

//./,j-.E„,.|-.",.».)}- 

J] (c + r{k, - ki)y^^^-'^^^ det I ^ I n dki. (E7) 
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Deforming the integral contour like kj — > kj + inj with Kj = LJj/2(5pj we obtain 



exp^E^^//•••/exp{-/3|:p,(fcJ-;.)}x 



\ 4/3p, 



n {c-in,-n,)+tik,-k)y^'^^'''^det ^ 



l[dki. (E8) 

ki^ki+iKi I 



It is important, that with the contour deformation no pole contribution is collected. (To 
see this we need the following observations. The first is, that in the different terms of 
det {dlj/dki}\^ _^^ ^^^ each (c — (kj — ki) + i{kj — ki)) has a power zero or —1, those terms 
with (c — {k,j — Ki) + i{kj — ki))~ cancel each other. The second is that if {kj — ki) > 
(thus c — {nj — Ki) can be negative), then Jjpi — Jipj > 1. The result is, that the con- 
tour deformation can be done such a way, that none of the zeros of the denominators are 
crossed.) Due to the prefactor in (]E8|) the order of magnitude of a general Ji, J2, ■ ■ ■ Jn 
term is exp | — X]j=i -Z^^</|/4/3pj [ times smaller, than that of the leading all J = term, 
i.e. supposing that pi is the largest element in the partition P, even the largest correction 
to flEll) means an 

relative error. As the possible maximum of pi is A^, this is well (over)estimated by the 
partition independent expression 

For those terms, which contribute to the partition function N/L is finite as A^, L — )■ 00, 
i.e. (lElOll decays exponentially with L, thus we may conclude, that the integral represen- 
tation of the sums over the quantum numbers is a very good approximation for the 6 Bose 
gas. 
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